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1 Introduction 


HERA data on Q 2 and x dependences of 9 ^q 2 , the logarithmic Q 2 derivative of the proton 
structure function F 2 (x, Q 2 ), have been published recently [I] The behaviour of d ^ A Q 2 
is of particular interest in the small x limit of deep inelastic scattering ( D1S ), where the 
DGLAP evolution equations || imply a relation 


dF 2 
d In Q 2 


2«s 
9 7r 


xG dglap (x,Q 2 ) . 


( 1 ) 


The most recent ZEUS data []2| are shown in Fig.l together with the corresponding 
GRV’94 H predictions. Each of these data points correspond to a different x and Q 2 
value. The ( x, Q 2 ) of Fig.l are averaged values obtained from each of the experimental 
data distribution bins. The wide spread of these ( x, Q 2 ) sets reflects the constraints 
imposed by the availability of very limited statistics in the limit of very small x. 

As can be seen, 9 c ^q 2 rises steeply with 4 up to x of approximately 5 x 1CT 3 , for 
Q 2 > 5 GeV 2 . This is in agreement with the prediction of perturbative QCD ( pQCD ) 
as manifested by the GRV’94 parameterization of F 2 (x,Q 2 ). The data shows a dramatic 
departure from this prediction in the exceedingly small x < 5 x 10“ 3 and low Q 2 < 

5 GeV 2 . In this kinematic domain, while GRV’94 predicts a continued increase of q^q 2 
with -, the actual data indicates that the logarithmic Q 2 slope of F 2 decreases significantly 
with -, up to the present experimental limit of x ~ 3 x 1CT 6 . 

There are some obvious experimental and theoretical observations and questions prompted 
by the above data: 

1. What is the detailed x and Q 2 dependence of d c ( A Q 2 ? As observed, the phenomenon 
of decreasing 9 ^q 2 is associated with small Q 2 and exceedingly small x. However, 
since the data is so sparse it is difficult to deduce its detailed behaviour . In partic¬ 
ular, the available data does not provide us with a clue as to the conduct of 9 ®*q 2 
at medium and high Q 2 , in the exceedingly small x domain. Note, that the entire 
data sample with x < 5 x 1CU 3 has Q 2 values smaller than 5 GeV 2 . 

2. The experimental behaviour of g^Q 2 as a function of x and Q 2 shown in Fig.l, 
indicates a possible departure from the DGLAP pQCD behaviour of Eq. i). for 
sufficiently small x. An interesting interpretation in this direction ( MRST ) has 
been recently proposed || assuming an initial “valence” gluon input distribution at 
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Figure 1: The HERA data and GRV’94 predictions (triangles) on the F 2 (x,Q 2 ) 
slope. The data are taken from Ref. |2|. 
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small Q 2 , which changes rapidly due to QCD evolution to the conventional “sea” 
distribution, at larger values of Q 2 . We shall comment on this below. 

The experimental F 2 (x,Q 2 ) data in the same kinematic region do not show any 
unconventional behaviour, i.e. F 2 is not sensitive enough, in the very small x domain, 
to resolve new features. The new data suggest that the logarithmic Q 2 slope of F 2 , 
is a more sensitive measure of possible deviations from the conventional pQCD 
dynamics. 


3. The unexpected phenomenon occurs in the transition region between “hard” and 
“soft” physics and one should carefully assess our ability to understand it within 
the framework of pQCD. Specifically, the predominant “soft” inclusive total real 
photoproduction cross section is experimentally finite. Namely, the Q 2 = 0 limit of 
a tot = 47r2 Q 2 em F 2 (x, Q 2 ) is finite at any W and thus F 2 oc Q 2 as Q 2 —■> 0. For a 
parameterization of the transition from DIS to real photoproduction, see Ref. [6]. 
However, the scale of this transition is not theoretically or experimentally specified. 
Generally, it is not clear if the observed change in the x and Q 2 dependences of 
q^q 2 can be analyzed utilizing pQCD techniques, or if this is just a signature for 
the dominance of “soft” non perturbative physics, when both Q 2 and x become 
small enough. Regardless, of its detailed dynamics, the data suggest a relatively 
smooth transition from the “hard” to the “soft” domains. Accordingly, we wish to 
match the “soft” ( non perturbative ) and the “hard” ( perturbative ) contributions 
probably around Q 2 ~ 1 — 2 GeV 2 , so as to better understand the “soft” limit of 
pQCD/ 


4. The theoretical interpretation of the new data presents a dual challenge . On one 
hand, it is desirable to understand and formulate the dynamics responsible for the 
changed behaviour of 9 ^q 2 , and to comprehend its implications for pQCD. On the 
other hand, regardless of this dynamics, the observed behaviour of 9 ®*q 2 implies 
that xG(x, Q 2 ) differs from our previous expectations in the limit of very small x 
and Q 2 . This is bound to impose changes on the input parton distributions which 
are currently used 0 , 0 , 0 - This is, indeed, the main point of Ref. || . 

5. An additional interesting study is to examine A = 9h Iq/iq from which one deter¬ 
mines the dependence of F 2 as a power of x at low x and given Q 2 , F- 2 cx x _A( R 2 ). 
A provides information pertinent to both Regge and pQCD analyses of F- 2 . As such 
we shall be interested in its Q 2 and x dependences, and limits, and the constraints 
imposed on A by unitarity [0. 


The goal of this paper is to expand and develop further the phenomenology and data 
analysis presented in our recent letter ||10|| in which the modified behaviour of 9 ^q 2 was 
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associated with the onset of unitarity screening corrections ( SC ).For an early suggestion 
that qJuQ 5 sensitive to SC , see Ref. [|lT] . 


To this end we study below the x and Q 2 behaviour of both 9 ^q 2 and check¬ 

ing the applicability of pQCD in x and Q 2 kinematic domains, not explored in detail 
previously. 


A reliable execution of this program depends on our ability to compare our calculations 
with the relevant data. As noted, such data is available only at a few isolated values 
of averaged x and Q 2 . This limited data base is not sufficient to test our suggested 
phenomenology. To overcome this difficulty, we propose to use the ALLM’97 computer 
generated data base which we call “pseudo data ”. The ALLM parameterization || 
for the proton DIS structure function is a 23 parameter description which provides an 
excellent reproduction of the a to t(l*p) data points. The updated fit 13 


points with 0 < Q 2 < 5000 GeV" and 3 x 10 


-6 


< x < 


to 1356 data 
0.85 yields a X 2 /ndf of 


0.97. Notably, ALLM’97 reproduces well both the recent data on 


dF 2 


and 


dlnF 2 


9 In 91n(l/a;) 

(this is shown later in Fig.|5|b). As such, we assume that ALLM’97 provides an accurate 
and reliable reproduction of the data, against which we can assess our calculations and 
determine our free parameters. 


The program of this paper is as follows: In Section II we briefly review the ALLM 
parameterization and its 1997 update. Section III is devoted to a theoretical review of the 
SC in hard parton physics and its relevance to F 2 , 9 {^q 2 and calculations. Our 

review examines the SC in the quark and gluon sectors. The systematics of 9 f*Q 2 and 

9hl(i/’r) anc ^ com Parison to the pseudo data, defined above, are discussed in Section 
IV. Our summary and conclusions are given in Section V. 


2 The ALLM’97 Pseudo Data Base: 


The ALLM parameterization 0 was introduced in 1991 so as to reproduce a to t (j'*p ), the 
total 7 *p cross section above the resonance region, for a very wide Q 2 range including 
real photoproduction ( Q 2 = 0 ). This is a multiparameter fit ( 23 parameters ) to all 
available data on a tot {xp) and F 2 , based on a Regge - type approach formulated so as to 
be compatible with pQCD and the DGLAP evolution equations. In its latest 1997 update 
a fit was performed to 1356 data points [|l], 0, for which an excellent x 2 fndf = 0.97 
was obtained. 

The ability to reproduce (T to t{l*p) over a wide kinematic range of 3 x 10 ~ 6 < x < 0.85 
and 0 < Q 2 < 5000 GeV 2 , is further demonstrated by the success of ALLM’97 in fitting 
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the F 2 logarithmic slopes, 9 ^q 2 and well. This is of particular interest, as the 


■ ■ x ) 

parameterization reproduces the F 2 and its logarithmic derivative data, also in the limit 
of exceedingly small x < 5 x 10~ 3 . 


The ALLM parameterization is presented by the following equations: 

Q 2 


F 2 (x,Q 2 ) = 


Q 2 + ml 


(f 2 f (x,Q 2 ) + F 2 R (x,Q 2 )) , 


( 2 ) 


where F R and F R are the respective contributions of the Pomeron and the secondary 
Reggeon exchanges to F 2 . 


F 2 p (x,Q 2 ) = C P (r)(^)^)- 1 (i _ x )PpC) 
F 2 R (x,Q 2 ) = 


(3) 

(4) 


ln Q +Q 0 

where r = ln | — C, ] is a slowly varying variable. 




The two modihed Bjorken x - variables are defined 


1 

xp 

1 

XR 


= 1 + 
= 1 + 


IV 2 —A / 2 . 
Q 2 + m'j, > 

W 2 -M 2 

Q 2 +m 2 


(5) 

( 6 ) 


M denotes the proton mass and mp and mp are htted effective Pomeron and Reggeon 
scales. The scale parameters ml, m 2 P , m 2 R and Ql control the smooth transition from 
Q 2 = 0 to high Q 2 . When Q 2 m 2 P and m 2 R , xp and xp approach x. Cp, ap, (3r and 

/ 3p increase with Q 2 as 

f(r) = fi + / 2 r /s (7) 

whereas Cp and ap decrease with Q 2 as 


9(r) = 9i + (gi ~ 92 ) [ 


1 + T 93 


!]• 


( 8 ) 


The parameterization depends on 23 parameters which are determined from the data 
fit. About half of these parameters are required to describe the low W ( high x ) region 
where higher QCD twist terms are important. A specification of the ALLM’97 fitted 
parameters is given in Table 2 of Ref. |T21. 


In general, Eqs. @ and 
duce two limiting cases: 


were constructed to be compatible with pQCD and repro- 
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(i) Reggeon - like behaviour at x —>• 0, with all corrections to the simple one Reggeon 
exchange absorbed in the dependence of the Reggeon intercept on r. 

(ii) The pQCD quark counting rules behaviour aX x —> 1, for which the power depen¬ 
dence of f3p and (3r on r are expected. 

In as much as the ALLM parameterization is constructed in a Regge-like scheme, 
it is important to note, that it is not a theory to be compared with the data, but a 
parameterization aimed at the best possible reproduction of the data. The ALLM’97 
parameters have two basic features dictated by the data, but not expected in a simple 
Regge-like theory. 

1) The Pomeron has a high scale of m 2 p ~ 50 GeV 2 . This high value is necessary to 
maintain a very smooth transition from predominantly “soft” Q 2 = 0 to a predominantly 
“hard” Q 2 of a few GeV 2 . A smaller value of m 2 P ~ 10 GeV 2 , suggested in the original 
ALLM fit [fj. implies a very fast transition, which is not compatible with the new data. 

2) Theoretically, for a single Reggeon exchange, a# should be independent of Q 2 . The 
Q 2 dependence, implied by the ALLM parameterization, secures the reproduction of the 
DIS data, indicating the need to correct the simple Regge picture in DIS. 

In this paper we evaluate the logarithmic derivatives g and derived from 

ALLM , Eq. (|2|) - Eq. The calculated 9 ®*q 2 and coincide with the experimental 

data in the kinematic region where they were measured. The ability to produce this 
detailed pseudo data base enables us to achieve two goals: 

1. To achieve a realistic reproduction of the behaviour of 9 c ^q 2 and over the 

complete Q 2 range, with a special emphasis on the small Q 2 and exceedingly small 
x domain. 

2. In our search for a dynamical explanation for possible deviations from the DGLAP 
pQCD expectations, we offer in section IV a pQCD calculation modified by SC which 
reproduces both the real data and the pseudo data points for Q 2 > 0.4 GeV 2 , as 
well as . 

A presentation and discussion of these issues will be given in section IV. 

An important item concerning ALLM is the uncertainty of the fit. An error calculation 
of a 23 parameter fit,where many of the parameters are highly correlated, is non trivial. 
Based on the estimates of Ref. || , we assess the pseudo data errors to be approximately 
8 - 10% . Over the kinematic region of interest the ALLM calculated errors are smaller 
than the comparable experimental errors. In section IV, which is devoted to the data 
analysis including the ALLM’97 pseudo data, we quote both errors whenever possible. 
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3 The Calculation of F 2 slopes with SC 


As stated in Ref.HTOl, we claim that a calculation of 


well for the data shown in Fig.l. 


a In Q 2 


including SC effects accounts 


I 11 the SC calculation we distinguish between two contributions leading to observed 
deviations from the DGLAP predictions of Eq. ([I]). These are: 


1. SC in the quark sector, corresponding to the passage of a qq pair through the target. 

2 . The analogous SC in the gluon sector, as xG(x, Q 2 ) appears in the expression for 
the slope of F 2 , Eq. 


3.1 Screening in the quark sector 


SC in the quark sector, calculated in the Eikonal approximation, were derived some time 
ago Jb|] |R| leading to an extensive phenomenology [[kj. In our context we have 


dF 2 sc (x,Q 2 ) 

dHQVQl) 

n Q {x,Q 2 ;b 2 t ) 


3 7 r 2 

27r«s 

w 


J db 2 {1 


_ e -KQ(x,Q 2 ;bt) 


}; 


r (b t ) xG 


DGLAP, 


x, Q 2 ) 


( 9 ) 

( 10 ) 


Ql is the photon virtuality scale from which we start the DGLAP evolution. Y(b t ) is the 
two gluon non-perturbative form factor of the target in impact parameter representation. 

1 


IW = 


; -i(bt-q t ) 


7T 


F(t) d 2 q t 


( 11 ) 


with t = —q 2 and F(t) denoting the form factor in t - representation. 

The fact that the b t - dependence factorizes from the x and Q 2 dependences was shown 
in Ref. [[TT| from the DGLAP evolution equations, using the factorization theorem |[TT|. 
It should be stressed that the factorization we use in Eq. ( |I0|) is only valid for t < Qq. 
Therefore, we assume below that only such small values of t are important in our integrals. 
This is the justification for simplifying the ^-dependence of Eq. (0), and using a Gaussian 
parameterization for Y(b t ): 

1 _ b t 


m) = 


R 2 


( 12 ) 
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Figure 2: The first order SC oc k 2 for F 2 (x,Q 2 ). 


where R 2 is a fitted parameter. 

In Fig.2 we present the lowest order SC to F 2 which are proportional to Kq, within 
the framework of the additive quark model ( AQM ) ( see Eq. (^) ). In the AQM we 
distinguish between two typical scales in the integration over q t . In the first diagram 
of Fig.2 we integrate over the distance between two constituent nucleon quarks. In the 
second diagram of Fig.2 the relevant integration is over the size of the constituent quark 
( see Refs. |L9| [^U| ). Eq. (^) assumes that the only constraint on the two gluons is their 
confinement within a nucleon whose size is R. 


The important property of Eq. (||) is that the calculated SC in the quark sector depends 
on a distance rj_ = -^ which we consider to be short enough for all values of Q 2 in our 

calculation. Consequently, we assess the perturbative calculation of the SC for 
to be reliable. This is not so for the calculation of the SC for F 2 , which has an important 
contribution from large distances r\ > ^2 as will be discussed below. As a result, the 
calculations of SC for F 2 are not that well defined, and may contain arbitrary errors. 


Following our previous publications [p~0|. plj], we define a damping factor 


dFj c (x,Q 2 ) 
p2 / \ _ 9In (Q 2 /<9o) 

Q {K) aF pGLAP (xQ2) ’ 

9 In Q 2 


(13) 


where ^(q^/q 2 ) calculated from Eq. The behaviour of Dq(k) as a function of x 
and Q 2 is plotted in Fig.3a , where we take K 2 = 10 GeV~ 2 . 
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3.2 Screening in the gluon sector 


The calculation of the SC in the gluon sector using an Eikonal approach was derived by 
Mueller |R| and further discussed in Ref.[^]. We obtain 

xG sc (x,Q 2 ) = ^ f Q d Q' 2 f dbti 1 ~ e- Ka{x '^' 2 ^ } + xG(x } Q 2 0 ) . (14) 

TC Z Jx X' JQl J 

Note that kg(x' , Q' 2 , b 2 ) = | k,q(x' , Q' 2 , b 2 ) defined in Eq. (P~0|) . 

We can now define a gluon damping factor 

xG sc (x,Q 2 ) 


D 2 g (x,Q 2 ) = 


xG DGLAP (x,Q 2 ) ’ 


(15) 


where xG sc (x, Q 2 ) is calculated from Eq. O. 

A difficulty in our calculation stems from the fact that the Q' 2 integration spans not 
only the short ( perturbative ), but also the long ( non-perturbative ) distances. Lacking 
solid theoretical estimates, we suppress the long distance contributions to Eq. (0) by 
imposing a cutoff on the Q' 2 integration, so as to neglect the contributions from rj_ > W 
integration. This ad hoc procedure makes our estimate of D G (x, Q 2 ) less reliable than 
the D‘q(x,Q 2 ) estimate. We have checked our cutoff procedure at Qq = 1 GeV 2 , using 
the GRV’94 parameterization for the solution of the DGLAP evolution equations and 
putting xG(x, Q 2 ) = 0 for Q 2 < Ql ( we shall denote this model GLMNO ) or imposing 
the GRV initial condition for xG(x,Q 2 ) at Q 2 = Qq ( GLMN1 ). The output result for 
Dq(x,Q 2 ) differs by about 10%, which lends credibility to our calculation. A graphical 
representation of Dq is provided in Fig.3b . Note that D 2 G is consistently smaller than 
Dq, i.e. SC in the gluon sector are bigger than the SC in the quark sector. We also 
observe that for any Q 2 the overall damping Dq ■ D G becomes significant once x is small 
enough. 


To check our sensitivity to the region of low Q' 2 we use a general property of the gluon 
structure function, namely, that gauge invariance requires that xG oc Q 2 at low Q 2 ( see 


Ref. |23|] for discussion ). Relying on this general property we assume that 


xG ( x , Q 2 < / l 2 ) 



(16) 


where /i 2 is a free parameter in the range 0.4 < fi 2 < 1 GeV 2 . We are motivated by 
the observation that even though GRV’94 evolves from a very low Qq = 0.4 GeV 2 , they 
actually fit the F 2 data only above Q 2 of about 1 GeV 2 . We have satisfied ourselves that 
our calculations are stable ( within 10 % ) to the choice of /j, 2 . 
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Figure 3: 


The damping factors: (a) Dq(x,Q 2 ), (b) Dq(x, Q 2 ) and 
(c) D 2 q (x,Q 2 )-D 2 g (x,Q 2 ). 
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3.3 Overall screening 


Our final expression for , i.e. the Q 2 logarithmic slope of F 2 with SC included 

in the calculation, is 


dFf c (x,Q 2 ) 

dHQVQl) 


D 2 q (x,Q 2 )D 2 g (x,Q 2 ) 


dF? GLAP (x,Q 2 ) 

dHQ 2 /Ql) 


(17) 


A detailed analysis of the systematics emerging from the above formula is given in the 
next section. 


The choice of input parton distributions obviously influences the output. At first sight, 
we have a very clear prescription of what to do. Indeed, the DGLAP evolution equations 
as well as our Eq. (|[J) and Eq. (fT3|) ( see also the full set of equations below, Eq. (p~8|) - 
Eq. (p2|), where the initial parton distributions are noted explicitly ) are written in such a 
way that we depend on the initial parton distributions at a fixed Q 2 = Ql, to solve them. 
These initial distributions, in principle, include everything that we have from “soft”, non- 
perturbative physics. All the information relevant to short distances rj_ < is included 
in the above equations both with or without SC. 


In an ideal situation, our procedure should be to use as input a set of parton distribu¬ 
tions which have been determined at Q 2 = Qq, without the influence of data at Q 2 > Qq, 
evolve to higher Q 2 using DGLAP evolution, and then correct for SC at Q 2 > Qy 
utilizing the damping factors. 

Unfortunately, such a set of input parton structure functions is not available, and we 
must make do with one of the available parameterization. These sets are obtained from 
global fits combined with DGLAP evolution equations, i.e. the determination of the input 
is also influenced by the data with Q 2 > Q q which may contain SC. We have chosen to 
work with GRV’94 parameterization since this is a relatively old parameterization and as 
such most of the fitted data has x > 1CU 3 . Our estimates show that in this kinematic 
region we expect the SC to be small. 


The GRV’94 parameterization has an additional advantage for us. It should be stressed 
that Eq. (H), Eq. (|L4]) and Eqs. (|T8|)-(|22|) in the following are proven in the double log 


approximation ( DLA ) in which we consider both as ln(l/x) > 1 and as \n(Q 2 /Q$) > 1. 
In the GRV’94 parameterization the initial value of Qq = 0.4 GeV 2 is so small, that the 
DLA is effectively valid above Q 2 ~ lGeV 2 . Indeed, above 1 GeV 2 GRV’94 provides a 
good reproduction of the DIS data. 


In a previous publication ]20| we have discussed the possibility that the proton is 
better described as a two radii object. This description is based on the large experimental 
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difference between the slopes in ^ of 7 p 
this option in Ref. 


J/'hpand 7 p —> J/'&X. We have formulated 


[20| and have applied it to gffXvQ 2 ) 
single and two radii models of the proton were examined. In as much as we endorse the 
two radii idea, we are unable, at this stage, to offer a decisive experimental signature that 
will discriminate between the one and two radii models ,when dealing with F 2 and its 
logarithmic slopes. This being the case, we limit ourselves in this paper to the simpler 
one radius approach. A detailed formalism for the two radii model was given in Refs. |20| 
and fid] 


in Ref. 


in which both a 


3.4 A screened calculation for A(Q 2 ) 


We list below the expressions for F^Qc^Q 2 ) and A (Q 2 ) = ^ which include 

screening corrections. The set of formulae are : 


Fi c (x,Q 2 ) 

dFj c (x,Q 2 ) 
d\n(l/x) 

Kq°(x, Q 2 ) 
xG sc (x, Q 2 ^j 

k» glap (x,Q 2 ) 


/?2 JU N f r Q2 

dQ' 2 

2vr- “ f Jqi 

{c + 111 k^ c (x,Q' 2 ) + E 1 (k% c (x,Q' 2 )) } 

+ Fr M (x,Ql) ; 

(18) 

732 AT N f r Ql 

dQ' 2 

2 7T 2 j± x f Jqi 

sinpc S0 M' 2 )) 1 

d\n(l/x) 1 J 

dFr ut (x,Ql) 

d\w{l/x) 

(19) 

Vv xGSCMi ) 

; (20) 

X f dx j c . 

7T Jx X JQq 

{ C + ln/cg GLAP (x / ,Q' 2 ) + E\ («g GLA V, Q' 2 )) } 

+ xG input (x,Q 2 0 ) ; 

(21) 

3 2 R2Q? XGDGLAF ^ Q2 ) ' (22) 


We have to integrate over the all distances including the long distance region. As we 
have discussed, this gives rise to undefined errors in our calculation which we attempt 
to estimate, assuming different cutoff values of the distances rj_ > A_. All integrals 
over the impact parameter b t , were done using the Gaussian parameterization for the 
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b t dependence of Eq. (0). In the above set of equations C denotes the Euler constant 
( C = 0.577 ) and Ei(x) = —Ei(—x ) is the exponential integral function (see Ref. |24||). 
N c and Nf are the colour and flavour degrees of freedom in QCD. Zf is the fraction of 
the electrical charge carried by a quark with flavor /. 


3.5 Asymptotic predictions 


It is instructive to discuss here the asymptotic predictions of Eq. (y) - Eq. (|22|) in the 
region of very small x. Such predictions have been discussed previously ( see Refs. 


10 


15, 19 


and we consider them here as a limit with which we can compare our actual 


calculations. 

In the region of small x, which is defined in our formulae as a region where 


ft, 


DGLAP 


(x,Q 2 ) = 


2vr a s (Q' 2 


xG 


DGLAP 


(x,Q 2 ) > 1 , 


VQ ' 3 R 2 Q 2 

the predictions for all observables which we discuss in this paper are as follows: 

1. Fi c (x,Q 2 ) =* ^ Ei f Z 2 Q 2 R 2 In k q (x,Q 2 ) ; 

? 2/02 


(23) 


2 . 


ahyQV’cp) ^ R- 2 Q 2 f° r x — x o, where Xq is determined from the equation: 
Ag GLAP (x 0 , Q 2 ^j = 1 . The solution to the equation Kq GLAP ( x, Q 2 ) = 1 is plotted 


DGLAP 


= 0.5 and 1.5. 


in Fig.4, where we plot also the solutions for k'q 
It should be stressed that at any value of Q 2 , for very small x < x 0 , the limit of 
the F 2 slope g^(Q^/Q' 2 ) independent of x, and is proportional to Q 2 . 

This limiting behaviour is reached at smaller and smaller x values as Q 2 is increased, 
and it may explain the Q 2 - dependence of the experimental data ( see Fig.l ) ; 


o dlnF^ c (x,Q 2 ) _ i 

’ <91n(l/ai) ^ ln(l/x) lnln(l/a;) 

dlnF^ c (x,Q 2 ) 


. Therefore, in the region of ultra small x: A (Q 2 ) = 
0, modulo logarithmic corrections independent of the value of Q 2 . 


d\n(l/x) 

We shall elaborate on this limiting behaviour in the next section ; 


4. xG sc (x,Q 2 ) ==> 2R2 ® 2 ln(l/a:) In k° glap (x, Q 2 ). Comparing this behaviour of 
the gluon structure function with the asymptotic behaviour of one can 

deduce that Eq. (HP does not hold in the region of exceedingly small x. 


Obviously, the above are just the exceedingly small x limits of the actual expressions, 
but they give us, together with the DGLAP predictions, a framework to discuss the 
experimental results, so as to develop a strategy of measurement which can test the 
asymptotic predictions. 
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Figure 4: Solutions to Kg = 0.5, 1 and 1.5. 
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4 Comparison with the experimental results 
and pseudo data 


The calculations presented in this paper were carried out with R 2 = 10 GeV 2 and 
Ql = 0.4 GeV 2 . We have not attempted to produce a “best fit”. However, we have 
satisfied ourselves that with these input parameters, we obtain good results which pro¬ 
vide a consistent reproduction of the measured and pseudo data. As stated, we have not 
pursued further the study of a two radii model jnj |2(J being unable to suggest a decisive 
signature which supports this hypothesis. 


The detailed features of our reconstruction of the ZEUS data [^| are shown in Fig.5a. 
As noted, the confusing compilation of experimental data points with widely spread 
(x, Q 2 ) values presented in this figure, makes it difficult to directly assess the x dependence 
of at a fixed Q 2 . Our analysis suggests that q^q^qI) grows monotonically with 

4 at fixed Q 2 , approaching its limiting value at x which is well below the x interval valid 
in our calculations. At high Q 2 and at fixed x, this rise is steep following the behaviour of 
xG DGLAP (x,Q 2 ) as expected from Eq. (!l[). In our approach, this is a direct consequence 
of the fact that in the high Q 2 limit and at fixed x, Dq Dq rs 1 throughout the x 
range of interest. Generally, for any fixed Q 2 we can find xq sufficiently small so that for 
x < Xq SC are important and, thus, both Dq and Dq are significantly smaller than 1 . 
However, from a practical point of view, as data are available for x > 2 x 10~ 6 and the 
GRV’94 parameterization is dehned for x > 1CT 5 , SC become relevant only at moderate 
and small Q 2 values ( see Fig.3 ). This general behaviour is reflected in our calculations 

where departs from the DGLAP predictions, and as Q 2 decreases the x depen¬ 

dence of g^Q^/Qiq becomes more moderate. Note that our approximations are not valid 
for x > 1CT 1 values. The final compilation of the calculated points, at (x, Q 2 ) values 
matching the experimental ones, reproduces the data very well, as is evident from the 
thick line of Fig.5a. The other lines in Fig.5a are for five typical constant Q 2 values, 
selected from the experimental data. Note that the thick curve in the figure, like the data 
points, belong to different Q 2 values. 


In Fig. 5b we show the same plot as in Fig.5a, only this time the fixed Q 2 behaviour 
is determined from ALLM’97 JT2|. Evidently, the behaviour of y, at fixed Q 2 , 

shown in Fig.5b is very close to our calculated behaviour presented in Fig.5a. 


Assuming the ALLM’97 pseudo data to be a reliable reproduction of the, unmeasured 
yet, real data - we can assess the various theoretical ideas and predictions for 9F q^q$ ^ , 
which are compared with the pseudo data in Figs 6 and 7. 
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9F 3 (x,Q 3 ) 9F 2 (x,Q 2 ) 

91n(Q 2 ) 91n(Q 2 ) 


(a) 





X 


Figure 5: (a) Our calculations for ZEUS data. The thick line is the prediction for 
the actual ZEUS data which were taken at values of Q 2 correlated with 
the value of x. (b) The same with ALLM’97 . 
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Figure 6: 


9F q\*q 2 ^ at fixed Q 2 . In addition to the ALLM band we show a typical 
data point with its error 
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A detailed comparison between ALLM’97 and the results obtained from our calculation 
( GLMN ) as well as GRV’94 || and MRST || are presented in Fig.6 for various Q 2 values. 
ALLM’97 is presented as a band which includes the estimated error We also show the 
typical small x experimental errors, so as to better assess the quality of the theoretical 
predictions. In Fig.7 we compare the pseudo data Q 2 dependence of 9F q^q 2 ^ at fixed x 
values of ICR 2 , 10 -3 , 10 -4 and 10 -5 with the same parameterizations, and the results of our 
calculation. Note that the MRST [|5] parameterization is valid only for Q 2 > 1.25 GeV 2 . 

1. Figs.5 - 7 show that the SC are able to account for the deviations from the DGLAP 
prediction, observed experimentally and reproduced in the ALLM’97 pseudo data. 
Our calculations are compatible with the observable scale Q 2 ~ 2 — 4 GeV 2 where 
such deviations start to be visible. 

2. In our approach we do not see any indication supporting an abrupt transition from 
the predominantly “soft” region to the hard one. This is compatible with the 
ALLM’97 fit which had to choose a high Pomeron mass scale m 2 p m 50 GeV 2 , 
so as to obtain a very gradual increase of the effective «p(0), with increasing Q 2 as 
required by the data. 

3. Figs.5 - 7 are compatible with a new scale of hardness ( Q 2 (x ) ) suggested in our 
main formulae ( Eq. ( fT8| ) - Eq. ( |22| ) ). Q 2 (x) is the solution of the equation 

k% glap (x,Q 2 (x)) = 1 . (24) 

The general features of this equation are shown in Fig.4, where the relevant Q(x) 
values can be identified. The novelty of this approach is that a hardness scale is 
derived for any x value under consideration. 

Comparing GLMN with GRV’94 and MRST we observed that both GLMN and MRST 
are significantly different from GRV’94 in accord with the real and pseudo data. We note 
also that in the exceedingly small x ~ 10~ 5 region the Q 2 dependence of MRST becomes 
different from ours but we are unable, at this stage, to check it experimentally. 

A summary of the Q 2 dependence of A (Q 2 ) = ^ at fixed x values of 10 -5 , 10~ 4 

and 5 x 1CV 3 is presented in Fig.8, where, once again, the ALLM’97 pseudo data is 
compared with the GRV’94 and MRST parameterizations and our results. The GLMN 
calculation was carried out with /r 2 = 0.6 GeV 2 controlling our continuation to very 
small Q 2 < 0.6 GeV 2 , where xG(x,Q 2 < fi 2 ) = ^-xG(x, // 2 ). To check our sensitiv¬ 
ity we plot in Fig.9 the results of several sets of calculations where in GLMNO we have 
xG(x,Q 2 < Qq) = 0 and Q 2 } — 1 GeV 2 , in GLMN1 we repeat the same calculation 
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putting the GRV initial condition for xG(x, Q%) at Ql — 1 GeV 2 . GLMN is as in Fig.8 
for three different values of /i 2 : fi 2 = 0.4, 0.6 and 1.0 GeV 2 . 

The study of A (Q 2 ) can be perceived from different points of view: 


1. In soft Regge approach we have, in the small x limit, oc x~ x where A = ap(0) — 
1 ~ 0.08. aep(O) is the “soft” Pomeron trajectory intercept at t = 0. In Regge 
approach the Pomeron intercept, and A, depend on Q 2 and x. Experimentally, the 
data available for 0 < Q 2 < 200 GeV 2 indicates a slow increase of A with Q 2 . This 
is also reproduced by ALLM’97 parameterization, where the increase persists up to 
the highest Q 2 ~ 10 4 GeV 2 described by this parameterization. This behaviour of 
A (Q 2 ) may be attributed to DIS dynamics or the existence of an additional “hard” 
Pomeron or both. 

2. The DGLAP evolution equations predict in the region of small x that 


A 


N 


7V e tt S lu(QVQg) 

ln(l/x) 


(25) 


i.e. we expect A to slowly increase with In Q 2 and to slowly decrease with ln(l/x). 

3. In the limit of exceedingly small x, SC predict that A oc independent of Q 2 . 

As can be seen in Fig.9, this limiting behaviour of A is reproduced in all the GLMN 
calculations where Eq. (|T6|) is used for very small Q 2 < /i 2 . 


Fig.8 suggests that for exceedingly small x, GRV’94 parameterization is consistent with 
ALLM’97 version, whereas GLMN and MRST are somewhat below and above ALLM’97 
predictions, respectively. Experimentally, there is no data to compare with. Moreover, 
it is difficult to assess if these are real differences or just numerical artifacts at the kine¬ 
matic edge of the various calculations. Note that both GLMN0 and GLMN1 agree with 
ALLM’97 for all x values, as shown in the figures. 

In our calculation for Figs. 8 and 9 one has to exercise extreme caution when evolving 
in 1 < Q 2 < 3 GeV 2 region, as the x dependence of is not yet asymptotic. 
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Figure 8: Comparison of different approaches: A (Q 2 ) versus Q 2 at fixed x . 
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Figure 9: Comparison of different models for the long distance behaviour of 
xG(x,Q 2 )\ A (Q 2 ) versus Q 2 at hxed x . 
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5 Summary and Conclusions 


In this paper we have attempted a detailed and systematic study of the proton structure 
function logarithmic derivatives ^ and Our calculations are carried out 

in the double log approximation of pQCD, including SC of the calculated quantities. The 
SC are calculated in the Eikonal model. Our study was motivated by the recent HERA 
results on 9F q^q 2 - showing a considerable departure from the DGLAP predictions in the 
kinematical domain of small Q 2 < 5 GeV 2 and x < 1CU 3 . 

A unique feature of our approach is that we compare our results with a pseudo data 
base, which is computer generated from the ALLM’97 parameterization. This is done so 
as to overcome the lack of detailed experimental data at sufficiently small x and Q 2 . This 
enables us to examine the fine details of our theoretical predictions, a task which can not 
be accomplished with the limited relevant experimental data presently available. 

The main conclusions of our study are: 

1. Our approach, in which we correct the unscreened DGLAP predictions for the effects 
of screening, enables us to achieve results on both logarithmic slopes of F 2 , which 
are in good agreement with, both the real and pseudo data, at x < 1CT 1 . 

2. The main feature of our investigation of 9F q^q 2 ^ is shown in Fig.[5| . It may appear 
that the particular structure of 9F g ^- , evident in the ZEUS data, indicates that 
once x is small enough then 9F q^q 2 ~ changes from an increasing to a decreasing 
function of K This is not the structure suggested by our calculations. We demon¬ 
strated that d ^)i^Q 2 ^ at fixed Q 2 both in our calculations, and in the pseudo data, 
remains a monotonic increasing function of -, in a good agreement with asymp¬ 
totic expectations. The SC only suppress the rate of such growth in comparison 
with the DGLAP approach. It is the combination of different (x, Q 2 ) data points 
which creates the particular structure seen in ZEUS data, and reproduced in our 
calculation. 

3. At low Q 2 we expect 9F q^q 2 ^ to be proportional to Q 2 , and therefore to decrease 
as Q 2 —> 0. This result is in full agreement with our calculations, as is evident from 
Fig-0 • Comparison with the MRST parameterization shows that this effect can 
be reproduced in the evolution equation, by suppressing the value of gluon (quark) 
structure function at low x in initial partons distributions. 
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4. Since the damping of d ~ significant only for small Q 2 < 5 GeV 2 and x < 
10” 3 , our calculations are on the boundary of being able to use pQCD. Moreover, 
this is a kinematic domain where both “soft” and “hard” dynamics are significant. 
Our calculations suggest a smooth transition between the “soft” and the “hard” 
domains, and as a consequence of this, the results that we obtain from “hard” pQCD 
calculations are shown to be stable, and compatible with the real and pseudo data. 
This observation is of a particular significance for the calculations of SC in the 
gluon sector, where our calculation also receives contributions from relatively long 
distances, for which pQCD provides no estimation of errors. We suggest that a scale 
of hardness be defined from the solution of Eq. (|23|). The above gives hope that 
the transition from ’’hard” to ’’soft” mechanism can be calculated theoretically. We 
consider this paper, together with our earlier Ref. [[TIJ , as a first attempt to quantify 
this description. 


5. Our study of A (Q 2 ) = dl ^n(i/x) ~ compatible with both the real and pseudo data. 
We confirm the expected asymptotic behaviour of A in the region of very small x 
(see Fig.[|). However, we would like to point out, that from Fig.|] one can see that 
the pseudo data indicate a different x dependence of A (Q 2 ) than our calculations. 
Only real data at low x can clarify the situation. 


6. Our general approach was to start from the unscreened DGLAP, for which we assess 
GRV’94 to be most suitable, and then correct for the effects of screening. An 
alternative approach would be evolving parton distributions, as done in ||. There 
are however, two main differences between our approach and MRST: the first is 
a practical deficiency of MRST, which is applicable only for Q 2 > 1.25 GeV" - 
this limit is somewhat high for the present analysis. The second is the different 
predictions between our approach and the MRST parameterization, in the region of 
small x ~ 1CT 5 (see Figs.|] - |j). 


We firmly believe that the experimental systematic studies of the F 2 slopes as well as 
other observables that are sensitive to the value of SC ( F^charm), Fl , — gi n q 2 ; etc- )> 
are needed to check one of the most fundamental problem of QCD: the value of the scale 
of the transition between perturbative and non-perturbative QCD. 
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